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1. INTRODUCTION 

In serious electric power systems market, there is a need to work electric utilities closer as far as 
possible. To analyze whether a power system can stay in a protected and solid working state under 
contingency conditions, those possibility cases that will bring about loss of load, loss of generation, or 
islanding are first to be distinguished. In view of framework administrators past experience each 
post-unexpected amount is allotted a level of seriousness as indicated by potential harm that could be forced 
on the power system by the amount, should the contingency happens [1]. 

In the current day power system planning and operation considerable interest is being shown in 
contingency analysis. Possibility screening and positioning is one of the significant segments of on-line 
system security appraisal. The target of contingency screening and positioning is to rapidly and precisely 
select a short rundown of basic contingency from an enormous rundown of potential contingency and rank 
them as per their seriousness. Reasonable preventive control activities can be actualized considering 
contingency that are probably going to influence the power system execution. Milano et al. [2] proposed an 
OPF-based market clearing calculation that represented voltage security limits. The security compelled 
optimum-path forest (OPF) program limits the pre-possibility target work while watching both the pre-and 
post-possibility framework imperatives [3]. A tale way to deal with multi target particle swarm optimization 
strategy for taking care of ideal power flow issue was proposed in [4]. Anitha et al. [5] an full depth 
reclamation (FDR) particle swarm optimization (PSO) technique proposed to take care of the OPF issue 
considering incline rate cutoff points of generators and line flow limits. Transformative calculations have end 
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up being ground-breaking meta-heuristic solvers when applied to complex issues [6]. Self-versatile models 
have demonstrated to improve the effectiveness of the calculations by and large. The self-versatile attributes 
are customarily given to calculations by controlling, under determination tension, the change rates utilized in 
the calculations. 

PSO was enlivened by the social conduct of feathered creatures running and was proposed to take 
care of persistent issues. The first PSO is a ceaseless PSO (CPSO) being applied and reached out to take care 
of numerous nonstop issues [7]. PSO has likewise been reached out to adapt to discrete issues [8]. Binary 
particle swarm optimization (BPSO) was initially evolved by Kennedy and Eberhart [9], to take care of 
numerous combinatorial issues, for instance, workshop booking [10] and feature selection [11]. In BPSO, the 
position is a twofold vector, and the speed is a nonstop vector. BPSO has been applied to some genuine 
issues. Sarath and Ravi [12] applied BPSO to produce affiliation rules from value-based datasets. The 
outcomes on a dataset from an Indian business bank demonstrated that applying BPSO gave more excellent 
principles as well as dodged repetitive standards. Taha and Nadi [13] utilized BPSO to identify the accessible 
frequencies in cognitive radio, which permitted to use the framework assets. 

The BPSO-based unique assignment accomplished higher recognition rates and lower bogus alert 
rates with various commotion proportions. Lin et al. [14] utilized BPSO to scan for exceptionally productive 
thing sets rather than visit thing sets in value-based databases. It was demonstrated that the BPSO-based 
calculation was increasingly proficient, progressively viable, and met quicker than GAs. In numerous 
different works, BPSO has likewise been changed to improve its presentation. In the first BPSO [15], a 
sigmoid function was utilized as an exchange work known as a shaped work. Mirjalili and Lewis [16], the 
V-formed transfer functions were proposed for BPSO. The position refreshing condition additionally thought 
about the past area. The trial results indicated that applying the V-molded functions while considering the 
past area improved the presentation of BPSO. Be that as it may, it was not satisfactory that which adjustment 
contributed more to the improvement. The exhibition of BPSO on various datasets intensely relied upon the 
particular transfer functions, in any event, when they were from a similar family (V-formed or S-molded). 
Consequently, it was anything but a simple undertaking to choose a proper transfer function for a specific 
errand or dataset. Islam et al. [17] proposed a period changing S-formed exchange work for BPSO which 
meant to give a smooth move from more investigation to more misuse during the developmental procedure. 
The trial results demonstrated that the time-changing transfer function helped BPSO to accomplish preferred 
outcomes over a standardized linear transfer function [18] and a notable V-molded transfer functions [16]. 
Liu et al. [19] gave a point by point investigation of the impact of the inertia weight parameter on the looking 
through capacity of BPSO. 

The BPSO algorithm is extremely important for practical and commercial use in the resolution of 
discrete problems, so this algorithm helps optimize power systems networks. The BPSO improves 
communication between populations and maintains population diversity in order to overcome limitations of 
classical optimization algorithms in resolving multi-parameter, powerful coupling and nonlinear problem 
optimization in engineering. The BPSO algorithm is based on the traditional PSO algorithm to balance global 
search and local search by adjusting the inertial weight. In this paper multiple objectives based optimum-path 
forest (OPF) method using binary particle swarm optimization technique was discussed. Fuzzy based 
contingency ranking method was used to find most severe lines of contingencies. The effectiveness of the 
BPSO method in giving solution for various objectives was tested on IEEE 30-bus and IEEE 14-bus systems. 


2. METHODOLOGY FOR BINARY PARTICLE SWARM OPTIMIZATION 

In BPSO, every arrangement in the populace is a binary string. Every binary string is of 
measurement n which is assessed to give parameter esteems. In the BPSO, every binary string speaks to a 
molecule. Strings are refreshed a tiny bit at a time dependent on its present worth, the estimation of that bit in 
the best (fitness) of that molecule to date, and the best estimation of that bit to date of its neighbours 
[20]-[25]. For binary strings, neighbors can be chosen in one of a few different ways. A few models are: (for 
an area of size k). Neighbors are the k paired strings whose Hamming separation is limited. For equivalent 
Hamming separations, the decisions are subjective. In the first place, self-assertively dole out gatherings of 
k strings to neighborhoods. Let the local size be the populace size. 

In standard PSO, everything is as far as a speed. By and large the speed is characterized as far as a 
likelihood of the bit evolving. In BPSO, a little bit at a time refresh is done probabilistically. At the end of the 
day, for a picked bit d in a picked string I it is changed to a 1 with a likelihood P that is an element of its 
inclination to be a 1, the best estimation of itself to date, and the best estimation of its neighbors. (1-P) is the 
likelihood of changing to a zero. When P is resolved, an irregular number R is produced. On the off chance 
that R<P, at that point the bit turns into a 1; else it turns into a zero. 
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Fundamentally, PSO and BPSO are the equivalent in the feeling of the consistent progression of the 
calculation, which in discovering pbest and gbest, update speed and position [8]. The main contrast among 
PSO and BPSO is the conditions used to characterize the updates of speed and position of every molecule. 
Equation (1) is utilized to refresh the speed while (2) and (3) are utilized to refresh the situation of every 
molecule. 


vig’) = wuiy + c11, (pbestia — Xia) + CoT2(gbestia — Xia) (1) 
sigmoid (vk) = —{ (2) 


1+e "id 


xk = F if rand < sigmoid (v%,) 3) 
iq 0, otherwise 

The x; and velocity vi of an i particle are described by its position. When searches are carried out, each 

particle modifies its location in line with its own experience and experience of all its surrounding particles, 

using the best position it has (pi) and everyone else (pg). 


3. PROBLEM FORMULATION 
Objective function-1: 


2 
F = DNS F; = DNS (a, PI + bP? + ci)Rs/hr 


where, F;is total cost function, p3 is generation of active power (control variables), NG is number of 
generator buses, ai, bi, ci: Coefficients of Fuel cost for i unit. 

Objective Function-2: Composite index based on line loadings, voltage profiles and voltage stability 
indices. 
Composite logic criteria (CLC)=(TIL_+TIve+TIvsi) 
CLC: Composite logic criteria 
TIu: Total index of line loadings 
TIvp: Total index of voltage profiles 
TIvs: Total index of voltage stability indices 

Figure | depicts the process of computing composite logic criteria. In this data obtained from 
Newton-Raphson load flow solution (line loadings, load bus voltage profiles and voltage stability indices) 
will be given to the program of fuzzy inference. The fuzzy system segregates the input data based on the 
rules framed for fuzzification and forms the total index for line loadings, voltage profiles and voltage stability 
indices. By adding these three indices composite logic criteria will be formed. As per objective function-2, 
composite logic criteria will be minimized, which in turn improves the security of the system. If the security 
of the system was improved then overloading in the lines will be reduced, voltage profiles of the load buses 
will be improved. Fuzzy rules to calculate CLC was given in Table 1. LL: loaded lightly, LU: loaded usually, 
CL: completely loaded, LO: loaded over, LV: low voltage, NV: normal voltage, OV: over voltage, VLI: very 
low index, LI: low index, MI: medium index, HI: high index, VHI: very high index, LS: low severe, US: 
under severe, SA: severe above, HS: higher severity, VLS: very low severe. 


Data of Line 
loadings , L= Th 
Iaa en Program of 
Data of Load Fuzzy m 
bus Voltages Inference > VP CLC 
Data of Voltage 
Stability Indices m~~ ~~~ Tivs 


Figure 1. Process of calculating composite logic criteria 
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Table 1. Fuzzy rule table 
Line Loadings Voltage Profiles Voltage stability Indices 
Input LL LU CL LO LV NV OV VU LI M H VM 
Output LS US SA HS US SA HS VLS LS US SA HS 


Above objective functions will be minimized subject to satisfaction of equality and inequality 
constraints. Newton—Raphson power flow methodology was used for power flow calculations. The equality 
constraint will be load balance equation i.e., total load and loss must be equal to total generation. The 
inequality constraints will be maintaining active power generation, reactive power generation, generator bus 
voltage profiles and load bus voltage profiles within their specified limits. Sets of non-linear flow of power 
equations, which regulate the system of power, are equality constraints. 


Fa Ete > IV, )¥,]e0s(@, - 6, +6,) = 0 (4) 
j=l 

Qoi -0n +Y NV, V, Y; sin( 0; -ô; + 6;) =0 (5) 
j=l 


In cases where, Pg; and Qgi are actual and reactive bus power outputs, load demand is represented by Pp; and 
Qp; at the same bus, and the bus admittance matrix components of |Y; jl and 6;; are represented. 

Unequal treatment constraints: These are the limitations representing the operating system and its 
security limits as follows: 
- Generators real and reactive power outputs 


pam < Pa SPRA eo NG (6) 
QG” < Qoi S QR i = 1.a, NG (7) 
- Voltage magnitudes at each bus in the network 
ymin < VV i = 1 NE (8) 
- Transformer tap settings 
pome pa Tr Sh aN (9) 
- Reactive power injections due to capacitor banks 
QURE < Qah, i= 1,..., CS (10) 
- Transmission lines loading 
§,< 57%, i =1,..,nl (11) 
- Voltage stability index 
Lj, < Li", i = 1,...,NL (12) 
4. FLOWCHART FOR IMPLEMENTATION OF BPSO 
Flowchart for the implementation of proposed methodology was shown in Figure 2. The main 
implementation steps as depicted in Figure 2 are given: i) initialize the multitude, the situation of particles is 
arbitrarily instated inside the hypercube of plausible space; ii) evaluate the exhibition of every molecule, 
utilizing its present position; iii) compare the presentation of every person to its best execution up until now; 


iv) compare the presentation of every molecule to the global best molecule; v) change the speed of the 
molecule; vi) move every molecule to another position; and vii) go to stage 2, and rehash until convergence. 
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Figure 2. Flow chart for implementation of proposed methodology 


5. RESULTS AND DISCUSSION 

Test cases in this paper are carried out on modified IEEE 30-bus and IEEE 14-bus systems. For 
IEEE 30-bus system the outage of line connected between buses 2 and 5 was obtained as top-1 severe 
contingency. For IEEE 14-bus system after conducting contingency analysis line connected between buses 2 
and 3 was observed as top-1 contingency. A comparison was made in between objective functions | and 2 
and the obtained results were reported in Table 2 for IEEE 30 bus system and in Table 3 for IEEE 14 bus 
system. From these results it can be observed that BPSO was effective in achieving specified objective 
functions. Figures 3 and 4 depicts comparison of active power generations under two objective functions for 
IEEE 30 bus system and IEEE 14 bus systems respectively. Figures 5 and 6 show comparison of Generator 
Bus Voltage Profiles under two objective functions for IEEE 30 bus system and IEEE 14 bus systems 
respectively. From these figures it can be observed that control variables are within limits. 


Table 2. Comparison of objective functions for IEEE 30 bus system 


Variable Objective Function 1 Objective Function 2 
Fuel cost ($/Hr) 828.4225 944.1078 
Composite index 1590 1347 
Loss (p.u.MW) 0.1392 0.0668 
Computation time (sec) 219.7030 371.61 
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Table 3. Comparison of objective functions for IEEE 14 bus system 


Variable Objective Function 1 Objective Function 2 
Fuel cost ($/Hr) 857.0573 1085.4 
Composite index 1016 957 
Loss (p.u.MW) 0.0867 0.0268 
Computation time (sec) 77.197 118.109 
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Figure 3. Comparison of active power generations under two objective functions for IEEE 30 bus system 


12- 
z — . 
= 14 —@— Objective Function-1 
ə 
a Objective Function-2 
c 
= 08 4 
= 
2 
= 
E o6 - 
c 
oO 
oo 
o 04 4 
> 
£ 
2 0.2 4 
y 
to 
Pl p2 P3 P4 P5 


Generator Number 


Figure 4. Comparison of active power generations under two objective functions for IEEE 14 bus system 
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Figure 5. Comparison of generator bus voltage profiles under two objective functions for IEEE 30 bus system 
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Figure 6. Comparison of generator bus voltage profiles under two objective functions for IEEE 14 bus system 


6. CONCLUSION 

From these results on the IEEE 30-bus and IEEE 14-bus systems, it can be seen that numerical 
results demonstrate the feasibility as well as the effectiveness of the OPF method. By generation rescheduling 
using proposed BPSO-OPF the security of the system can be improved effectively. This BPSO-OPF method 
is useful for providing an efficient solution for multiple objective functions. BPSO based OPF will be a very 
useful tool to operate and manage the electric power systems in efficient ways and it can be used for 
optimization in practical power systems. The results clearly indicates that the OPF method works 
satisfactorily in bringing a superior solution for the specified objective function while maintaining all control 
variables within limits. 
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